Pesin-Type Identity for Weak Chaos 
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Pesin's identity provides a profound connection between entropy Iik s (statistical mechanics) and 
the Lyapunov exponent A (chaos theory). It is well known that many systems exhibit sub-exponential 
separation of nearby trajectories and then A = 0. In many cases such systems are non-ergodic and 
do not obey usual statistical mechanics. Here we investigate the non-ergodic phase of the Pomeau- 
Manneville map where separation of nearby trajectories follows 5x t = Sxoe Xat with < a < 1. 
The limit distribution of A Q is the inverse Levy function. The average (X a ) is related to the infinite 
invariant density, and most importantly to entropy. Our work gives a generalized Pesin's identity 
valid for systems with an infinite invariant density. 
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Chaotic systems are characterized by exponential sep- 
aration of nearby trajectories, which is quantified by a 
positive Lyapunov exponent A [H . Such a behavior leads 
to the need for statistical approaches since chaos implies 
our inability to predict the long time limit of a system 
in a deterministic fashion. A profound relation between 
chaos and statistical mechanics is given by Pesin's iden- 
tity [l[ . It states that the Kolmogorov-Sinai entropy h KS 
is equal to the Lyapunov exponent A for a closed ergodic 
Id systems (to the sum of positive Lyapunov exponents 
for d > 1). At the same time, it is well known that many 
systems such as Hamiltonian models with mixed phase 
space S, systems with long range forces certain bil- 
liards P| and one-dimensional hard-particle gas [B| have 
a Lyapunov exponent equal zero. While for complex sys- 
tems it may be extremely difficult to determine whether 
the Lyapunov exponent is zero or small, due to numeri- 
cal inaccuracies, it turns out that most fundamental text 
book examples of chaos theory may have a zero Lyapunov 
exponent. Prominent examples for such weakly chaotic 
systems are the logistic map at the edge of chaos (Feigen- 
baum's point) Q and the Pomeau-Manneville map which 
is used to model intermittency (originally in turbulence) 
Q. If the Lyapunov exponent is zero, i.e. separation 
of trajectories is sub-exponential, we have a strong in- 
dication that the usual Boltzmann-Gibbs statistical me- 
chanics is not valid. Indeed it was found that certain 
systems with zero Lyapunov exponents break ergodicity 
d, [§] . Classical entropy theory is also not applicable in 
this case [f| 0] , particularly the entropy and average al- 
gorithmic complexity grow non linearly in time 0] , while 
for a system with a positive Lyapunov exponent they in- 
crease linearly in time. Still the situation is not hopeless 
from the point of view of statistical mechanics and one 
may consider distributions of time average observables 



Connection between possible generalizations of usual 
statistical mechanics and weakly chaotic systems have 
attracted much attention recently. In particular, a gen- 
eralized Pesin's identity for the logistic map at the edge 



of chaos was investigated using Tsallis statistics [g, 1 131 . 
A critical discussion of this approach is given in Ref. jl4| 
(and see a reply in [HI). According to [3] a meaningful 
generalized Pesin's identity must satisfy certain require- 
ments, (i) Averages must be made with respect to the 
natural density, in our case the infinite invariant density 
3 12 1 (see details below), (ii) Evolution should be for 
long times, unlike previous attempts to generalize Pesin's 
identity, (iii) The relevant entropy is the entropy of Kol- 
mogorov and Sinai, not Boltzmann-Gibbs [18|. (iv) Re- 
sults should be general in that they do not depend on par- 
ticular initial conditions. The generalized Pesin's iden- 
tity, Eq. (frnj) below, fulfills these requirements. Thus, 
we establish a profound relation between separation of 
nearby trajectories and entropy, even though the separa- 
tion is sub-exponential. 

Consider the Pomeau-Manneville map [iji] on the unit 
interval with one marginally unstable fixed point located 
at x = 

M(x t ) = x t + ax\ (mod. 1) , z > 1, a > 0. (1) 

The discontinuity point £ is defined by M(£) = 1. This 
map is one of the pioneer models of intermittency. Its 
generalizations attracted vast research using different 
methods such as continuous time random walks 2(| 21 1 
and periodic orbit theory [22l | to name a few. Sojourn 
times of trajectories in the vicinity of the unstable fixed 
point are described by power law statistics leading to 
aging and non Gaussian fluctuations 0|, which are 
related to weak ergodicity breaking [H, @ . 

In the non ergodic phase z > 2 the density function 
of the map is concentrated on the unstable fixed point 
in the long time limit. The derivative |M'(x)| at this 
point is equal to 1, so A = as shown already in 0]. 
Such a behavior is found since most of the time the par- 
ticle spends in the vicinity of the marginally stable fixed 
point. Following [7| assume that the sensitivity of nearby 
trajectories is stretched exponential Sx t = fao^"'" with 
< a < 1. Using the chain rule, and the dynamical 
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mapping x t +i = M(x t ) we have 



1 t_1 

Aa(^o) = -^ln|M'(x i )| 



(2) 



where the dependence on initial condition is emphasized. 
For normal case and an ergodic system we have a = 1. 
Then the usual Lyapunov exponent is [l[ 

t-l 

A= lim - Vln|M'(xi)| = / da; p(:r) In \M'(x)\ , (3) 

where p(x) is the invariant density of the system. Er- 
godicity ensures that the time average is equal to the 
ensemble average. To prove that a system actually ex- 
hibits stretched exponential separation of trajectories, it 
is sufficient to find the limit distribution of A a (iro) and 
show that it is not trivial (i.e. < a < 1). As we now 
show in the non ergodic phase X a does not converge to a 
constant but remains a random variable. Below we ob- 
tain the distribution of X a and for this we now calculate 
the density of trajectories p{x,t) of Eq. (TT]). 

To obtain the density analytically we use the approx- 
imation of the map [24j for i < 1, dxt/dt ~ ax*, and 
extend it to be valid on the interval (0, £). When the 
trajectory reaches the boundary x = £ it is randomly 
reinjected back to the interval (0, £). The density func- 
tion p c (x,t) of this system is governed by the equation 

= _J_ (ax*p c (x,t)) +at*p c (Z,t), (4) 

where the subscript c in p c {x, t) is for continuous approxi- 
mation. The first term on the LHS represents determinis- 
tic escape from the marginally unstable fixed point while 
the second term accounts for reinjection of particles. The 
solution of Eq. (TJJ in Laplace space is 



Pc(x,s) 

where 

6 x (s)=b(z-l) 



(5) 



(6) 



and b = (z — l)~ 1 a~ 1 x 1 ~ z . Statistics of the system is 
controlled by a = 1 for z < 2, a = — ^j- for z > 2. 
Considering the small s behavior (equivalent to t — > oo) 
and transforming the solution into the time domain, we 
obtain for < a < 1 (z > 2) 



p c (x,t) 



i(7ra) 



t a - 



sin(7ra) , a 

For z < 2 we find the following solution 



x ^ 
x = 0. 



Pc(x,t) 




x^O 
x = 0, 



(7) 



(8) 




FIG. 1: Numerically calculated p(x) — t 1 ~ a p(x, t) for the map 
Eq. |T} with a = 0.3 and a = 1. Here ti = 10*, from bottom 
to top i = 2,3, 4. Dashed lines correspond to Eq. (|10|l with 
no fitting parameters. x c represents the crossover from one 
asymptotic to another. Note, that x c decreases with time, 
and when t — > oo we approach the infinite invariant density. 



and for z = 2 the solution is given by 



^0 
a; = 0. 



(9) 



Note, that the density function is time independent only 
for z < 2 and x ^ Eq. J5J. We introduce the infinite 
invariant density 



1 sin(Tra) 



p c {x) =t'- a Pc (x,t) 



i(7ra) 



x^Q 
x = 0, 



(10) 

for < a < 1. Scaled functions /O c (x) are independent 
of time for x ^ 0. Note, that Pc^) ~ x _1 / Q , and its 
integral diverges, dxp c (x) = oo. Thus, p c (x) is not 
normalizable [lg, [l7|- Still, as we show later, the infi- 
nite invariant density is useful for the calculation of the 
statistical properties of the dynamics. 

We compute the invariant density numerically and 
compare it with the analytical density function Eq. (fTU)) . 
In these simulations we start with a uniform density and 
plot p(x) — t 1_Q p(x, t) versus x. Results are shown in 
Fig. Q] We find excellent agreement between Eq. (fl0|) 
and numerics without fitting. Horizontal lines represent 
asymptotic solution for x = calculated for the corre- 
sponding time of the simulation, while the sloping line 
corresponds to the asymptotic solution for x ^ which 
decays as x" (see Eq. ([T0|) ). In Fig.^jxc represents the 
crossover from one asymptotic of p{x) to another. As 
t — > oo, x c — ► and we approach the infinite invariant 
density. 
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To find the distribution of X a (xo) Eq. |2]), we use a sim- 
ple stochastic approach. The same distribution can be 
found using Darling-Kac theorem applied to our observ- 
able Eq. [T2|. Consider the logarithm of the deriva- 
tive of the map y — In \ M'{xt)\ in the non ergodic phase 
< a < 1. We define a two state process I(t) = if 
y < £ and I(t) = 1 if y > £. Waiting times in state are 
distributed according to the probability density function 
(PDF) tp(t) ~ A/t 1+a as t — > oo, or in Laplace space 
0(a) = J °° dte- st ip(t) ~ 1 - Bs", as s -> 0, where A, B 
are positive constants, so the average waiting time is in- 
finite as well known [2C 



21| . In contrast, waiting times 



in state 1 have a characteristic average time. Neglecting 
correlations, we consider I{t) as a renewal process. Let 
n be the number of renewals, namely, number of transi- 
tions from state to 1. The logarithm of the derivative 
of the map In |M'(cc t )| is equal to zero most of the time 
(roughly when I(t) = 0) since the trajectory stays for 
long time near marginally unstable fixed points, only for 
short periods its value deviates from zero. The sum of 
logarithms along a trajectory is thus proportional to n: 
Si=o m \M'(xi)\ ~ c n, where c is a positive constant. 
Our goal is to calculate the PDF of scaled generalized 
Lyapunov exponents A Q Eq. (|2|) . The PDF of the num- 
ber of renewals n which occur up to time t is given by 

M 



(Bn) 



l/a 



(11) 



where l a is the one-sided Levy PDF defined through 
its Laplace transform l a (s) = exp(—s a ). We define 
C = A a / (A Q ) and, since (X a ) — c (n), ( = n/ (n) is inde- 
pendent of c. Using Eq. (JTTJ) 



p*(0 



r l/a(! 



aC, 1+1/c 



In 



F 1 / a (l + Q;) 
(5> 



(12) 



This is one of the main equations in the manuscript since 
it gives the PDF of scaled generalized Lyapunov expo- 
nents A a /(A a ). Distributions of £ = A Q / (A a ) obtained 
by simulations are shown in Fig. [21 Smooth curves cor- 
respond to analytical PDF Eq. (fT2"]) without fitting. The 
perfect agreement between theory and numerical results 
indicates that the general theory works well for finite time 
simulations. 

Now we calculate the average (A Q ). Using Eq. ([2]) 



(Ac 



p(x ) dx 0l (13) 



where the averaging is over initial conditions distributed 
according to some initial density. Since we are interested 
in the long time limit, we replace the summation with an 
integral and average over the density function 

(A Q )~i/ dx f \n\M\x)\ p{x,t)dt. (14) 
1 Jo Jo 



Z£ 0.4 



Q 
0. 




FIG. 2: P a (() for a = 0.75,0.59,0.5 from bottom to top on 
the LHS of the figure obtained with (A a ) calculated according 
to Eq. (|16p . Here t — 10 5 . Smooth curves correspond to 
analytical PDF Eq. Jl2j without fitting. 



According to Eq. Q, the density function has two 
asymptotics valid for x = and x ^= 

< A «> ~ 7S I ' dx [ p(0,t) dt+ 

t Jo Jo 



+ — dx \n\M'(x)\ p(x,t) dt, (15) 



1 

where x c denotes the crossover from one asymptotic of 
p(x,t) to another (see Fig. [T|). We define it as p(x c ) = 
p(x = 0). Using Eq. 0, x c = a a t~ a . When t -> oo, 
x c — > 0, so the first integral in Eq. (| 1 5j) vanishes. This 
fact is not obvious since p(x = 0) — > oo (see Eq. (0), it 
happens because we consider a specific observable with 
In |M'(x)| — » as x — > 0, which cancel the t a divergence 
found in Eq. I|7|). Computing the integral over time and 
using t p(x, t) = p{x), we obtain our second main re- 
sult 



1 



dxln|M'(a;)| p(x). 



(16) 



Thus, even though p(x) is not normalizable it yields 
the average generalized Lyapunov exponent (A Q ). Since 
ln|M'(x)| vanishes precisely where the infinite invariant 
density diverges, the integral is finite and positive. Our 
main result Eq. (| 16[) is very elegant, since up to a con- 
stant a it states that all one needs to describe separation 
of trajectories is to replace the invariant density with the 
infinite invariant density. For the stochastic model Eq. 
with < a < 1 using Eq. ([TO]), we find 



(X a ) - - 



dx 



sin(7ra) ln(l + azx<* ) 



-l/a 



(17) 
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z 



Our approach opens up many possibilities of finding 
connections between dynamical systems characterized by 
weak chaos and weak ergodicity breaking. Investigation 
of the relation between weak chaos and transport, gen- 
eralizing the escape rate formalism [l[ is one example. 
Similar to the generalized Lyapunov exponent X a found 
here, we expect that time averaged transport and diffu- 
sion coefficients remain random variables as was recently 
shown in stochastic models of weak ergodicity breaking 
[2^ . Randomness of A Q indicates the modification of 
the whole structure of statistical mechanics for weakly 
chaotic systems. 

This work was supported by the Israel Science Foun- 
dation. We thank R. Klages, P. Howard, A. Robledo and 
R. Zweimuller for discussions. 



FIG. 3: Numerical average of (X a ) (triangles up) perfectly 
agree with (X a ) (triangles down) calculated by Eq. (I16|l when 
p is calculated numerically. Here t = 10 5 . Solid line is (X a ) 
Eq. (|17[) found with approximate infinite invariant density 
function. 



We emphasize that our main results Eqs. (jT21|T5|) are gen- 
erally valid for systems with an infinite invariant density 
p(x) . Eq. (fT7|) is specific to maps with one unstable fixed 
point and can be used to verify the theory numerically. 

In Fig. [3] we present perfect agreement between numer- 
ical simulations of (A Q ), and Eq. (fl"6|) with p{x) calculated 
numerically. For not too large z, good agreement between 
simulations and the theory based on the stochastic ap- 
proximation for the infinite invariant density Eq. (|I7I) is 
found. For large z the convergence is slowed down, since 
a is small. In the ergodic phase z < 2 the standard 
Lyapunov exponent Eq. ([3]) is recovered. 

For weakly chaotic systems, (A Q ) can be used to extend 
Pesin's identity. Mathematicians have rigorously shown 
that entropy for the maps with infinite invariant measure 
h a satisfies Rohlin's formula 17, 2(| 27, 28| 



dx p{x) In |M'(ie) 



(18) 



where p{x) is the infinite invariant density. From Eqs. 
(fTBl [T8]) we obtain the identity 



h a = a (A Q ) 



Thus, we have established a profound relation between 
statistical mechanics and chaos theory for weakly chaotic 
systems in terms of Pesin-type identity Eq. (fl9)) . The 
connection between (A Q ) and h a is in fact deeper than Eq. 
(fT9|) . because according to Eq. (TT2"| the distribution of the 
generalized Lyapunov exponents A Q is controlled by the 
scale h a since C = X a / (A Q ) = a\ a /h a . Only in the limit 
a — > 1 from Eq. (Ti~2"j) we get lim a _vi P a (X) = 6(X — hxs)- 
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